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We study a crossover from classical to quantum picture in the electron energy statistics in a 
system with broken time-reversal symmetry. The perturbative and nonperturbative parts of the 
two level correlation function, R(iv) are analyzed. We find that in the intermediate region, A <C 
w ~ t^ 1 <C tJrg, where tE and t erg are the Ehrenfest and ergodic times, respectively, R(uj) consists 
of a series of oscillations with the periods depending on tE, deviating from the universal Wigner- 
Dyson statistics. These Ehrenfest oscillations have the period dependence as t^ 1 in the perturbative 
part. [For systems with time-reversal symmetry, this oscillation in the perturbative part of R(u>) 
was studied in an earlier work (I. L. Aleiner and A. I. Larkin, Phys. Rev. E 55, R1243 (1997))]. 
In the nonperturbative part they have the period dependence as (A -1 + crfe) with a a universal 
numerical factor. The amplitude of the leading order Ehrenfest oscillation in the nonperturbative 
part is larger than that of the perturbative part. 
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I. INTRODUCTION 

In the last decade, there has been increasing interest 
in investigating interference phenomena of electronic mo- 
tion in ballistic quantum dots 1 . In semiconductor quan- 
tum dots the confining potential has the length scale 
much larger than the Fermi wavelength, and therefore 
one may expect that electron motion in such systems may 
be described by methods of classical trajectories. Typi- 
cally, the classical electron motion inside dots is chaotic. 
An important question is, what happens to the energy 
spectrum. It turns out that the energy spectrum of a sin- 
gle electron in chaotic quantum dots is highly sensitive to 
the parameters of the system (e.g., Fermi momentum, the 
strength of the applied magnetic field, the impurity con- 
figuration, the shape of the boundary of ballistic quan- 
tum dots, etc.). Substantial progress has been made by 
studying the statistics of energy levels. The level statis- 
tics has played an important role in the theories of atomic 
nuclei 2,3 , disordered metals 4-6 , and quantum chaos 7 . A 
particularly interesting quantity in the level statistics is 
the so-called two level correlation function, R(u>): 

R(u) = A 2 (Tr(e + | - H)Tr(e - | - H)) (1) 

where A is the mean level spacing and H is the Hamilto- 
nian. (■ ■ •) denotes the average over a wide energy band 
(for individual quantum dots) or the impurity configura- 
tion (for quantum disorders). 

In his pioneering work, Dyson divided the Hamiltonian 
matrices into three classes: Gaussian orthogonal ensem- 
ble (GOE), Gaussian unitary ensemble (GUE) and Gaus- 
sian sympletic ensemble (GSE) and used the random ma- 
trix theory (RMT) to find that R(u) follows different 
universal behavior for these ensembles 8 ' 9 . Gor'kov and 
Eliashberg were the first to use Wigner-Dyson ensemble 
to study transport properties in small metallic grains 4 . 



For disordered metals, the behavior of R(oj) is now well 
understood. For grains in the presence of a lot of quan- 
tum (Born) impurities, theoretical justification of the 
universal RMT result came up with the supersymmct- 
ric field theory due to Efetov 5 , assuming that u> is much 
smaller than the Thouless energy Erh = hD/L 2 with D 
being the diffusion constant and L the size of the grain. 
The behavior of R(u>) beyond Exn is nonuniversal 6 . That 
is, the level statistics at w ^ Exh is system-dependent. 
In Ref. 10, it was found that there is a small constant 
correction to R(u>) at u> ~ A due to the nonuniversal be- 
havior of electron motion. Recently, the weak localization 
correction to R(uS) beyond the Thouless energy Exh was 
studied in Ref. 11 with the help of replica technique 12 . 

Bohigas, Giagonni, and Schmit (BGS) proposed a con- 
jecture that the fluctuations of levels of quantum chaotic 
systems may follow the universal RMT results 13 . One 
routine method to study the behavior of R(u>) in ballis- 
tic quantum dots is to employ the Gutzwiller formula 14 . 
A cornerstone was put by Berry 15 . By taking into ac- 
count the diagonal contributions, he was able to use the 
Hannay-Ozorio de Almeida (HOA) sum rule 16 to repro- 
duce the leading perturbative term of the universal RMT 
results. Proceeding along this line, an important progress 
was made in Ref. 17, where the off-diagonal contributions 
were studied and the oscillation similar to the universal 
RMT results was thereby found. 

In this paper, we will study the ballistic quantum dot 
with weak diffractions, where the classical electronic mo- 
tion is chaotic (chaotic quantum dots). In such case, 
there is a new scale 4e, so-called Ehrenfest time, which 
has logarithmic dependence on h: 

t , = I|m(£)|~i|lnA|. (2) 

Here A is the Lyapunov exponent, L is some macro- 
scopic size and Xf is the Fermi wavelength. Physically, 
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it is the time scale for an initial Gaussian wave packet 
(with the typical spatial spreading ~ Xp) to develop to 
some macroscopic size. At this scale, the electron mo- 
tion shows a crossover from classical to quantum pic- 
ture. Such a crossover exists in other phenomena also. 
Larkin and Ovchinikov 18 found that it is possible to use 
the method of classical trajectories in the theory of su- 
perconductor for t < tE only. Besides, they were able to 
estimate tE for Lorentz gases of semiclassical impurities. 
The logarithmic dependence on h of t E was also found 
by Berman and Zaslavsky in the studies of classical chaos 
of the kicked rotator model 19 . In the past two decades, 
it has been found that the Ehrcnfcst time plays an im- 
portant role in condensed matter. Wilkinson addressed 
the importance of such a scale in the semiclassical stud- 
ies of sum rules over generic matrix elements 20 . Arga- 
man studied the conductance in the scale lj -C t E 121 . In 
Ref. 22, Aleiner and Larkin not only gave the estimation 
for t E , but also studied the weak localization correction 



to the conductivity at the crossover scale lj - ~ o E 
showed that diffusons do not couple at time t <t E while 
at time t>t E , the coupling reaches a universal value. As 
a result, the amplitude of the coupling between diffusons 
(Hikami box) 23 is an oscillating function of Lots- Later 
it was found that such amplitude has important applica- 
tions in studying the crossover from classical and quan- 
tum picture of shot noise 24,25 and the density of states 
(DOS) in Andreev billiard 26 - 27 . 

For chaotic quantum dots with one macroscopic scale 
L, the ergodic time t erg ~ L/vf, which is of the same 
order of A -1 . In the semiclassics, t E opens a new region: 
A <C t E x ~ lj <C t~}g- The manifestation of tE in the 
level statistics of chaotic quantum dots was first stud- 
ied in Ref. 28. According to the RMT, asymptotically, 
for large lj the two level correlation function consists of 
two parts: One part is perturbative in — . The other is 
nonperturbative, which oscillates with the period A. In 
Ref. 28, it was found that there is a small oscillation with 
the period proportional to t E x correcting the perturba- 
tive part of the universal R(ui) for GOE 29 : 



A 3 d 2 



-ilj 



where T 2 (oj) is 



T 2 (lj) = cxp iLots 



(3) 



(4) 



Here A is the Lyapunov exponent and A2 characterizes 
the fluctuations of A [see Eqs. (A15) and (A16) for the 
general definitions] . Typically, A and A2 are of the same 
order. 

In this paper, we study the two level correlation func- 
tion in the system with broken time-reversal symmetry. 
We recall that in this case, the RMT predicts the famous 
result that the two level correlation function is described 
by the elegant formula: 



R u (lj) 



A 2 



2tt 2 c 



1 — cos ■ 



(5) 



Asymptotically, for large lj, the two level correlation 
function is exactly truncated at the (A/lj) 2 term 5 ' 8 . 
An important question is whether it is the property of 
generic systems with broken time-reversal symmetry or 
just comes out with the RAIT for GUE. In this paper, 
we will show that Eq. (5) just comes out with RMT, or 
describes the universal limit of the two level correlation 
function in quantum disorders. For ballistic quantum 
dots, surprisingly, we find that in the crossover region 
t E x ~ lj <C t~rg, the behavior of the two level correla- 
tion function deviates from the universal Wigner-Dyson 
statistics given by Eq. (5) . More exactly, such deviations 
are formulated as 



AR»(lj) = AiZ» + AiC (lj), 



They where 



Ai?» = -—Re 
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T 3 (lj) = cxp ( -lLjt E - ^-^2 
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(9) 



Compared to Eq. (3), the deviations in systems with 
broken time-reversal symmetry are more interesting. 
First, we point out that the main feature of these de- 
viations is the appearance of two types of oscillations 
in addition to the oscillation described by the Wigner- 
Dyson statistics. One type correcting the perturbative 
terms has the period ~ t^ 1 , which has the similar origin 
as Eq. (3). In fact, as we will show in Sec. Ill, Eq. (3) 
comes from the two loop approximation. In the bro- 
ken time-reversal system, the two loop approximation 
vanishes. Instead, Eq. (7) comes from the three loop 
approximation. They arc the weak localization correc- 
tions to leading universal R(lj) in the crossover region. 
Eq. (9) belongs to the other type of oscillation, which 
corrects the nonperturbative part of the universal R(lj) 
described by the universal Wigner-Dyson statistics and 
was not found in the previous work 28 . It has the period 
(A -1 + at E ) ■ Both of them have ^-dependent peri- 
ods. In this way, we term them the Ehrcnfcst oscillations. 
In the broken time-reversal system, the oscillation with 
the period (A -1 + at E ) is stronger than the oscilla- 
tion with period ~ t E x . For the time- reversal system, 
the Ehrenfest oscillation with the period (A -1 + at E ) 1 
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does exist, but the amplitude is smaller than that de- 
scribed by the Wigner-Dyson statistics, which is even 
small. 

In Ref. 30, it is proposed that the amplitudes of both 
perturbative and nonperturbative parts of R(u) are re- 
lated to each other through the classical spectrum deter- 
minant, which depends on the eigenvalues of the Perron- 
Frobenius operator. In the universal limit lo <C t~^ g , it 
coincides with the result of RMT. Indeed, for the broken 
time-reversal system, we see from Eq. (5) that the ampli- 
tudes of both parts are in the same order. However, such 
a conjecture docs not exist once the quantum corrections 
are taken into account. In fact, the leading quantum cor- 
rection in the perturbative part, i.e., Eq. (7) comes from 
the three loop approximation, while the leading quan- 
tum correction in the nonperturbative part, i.e., Eq. (9) 
comes from the two-loop approximation. As a result, the 
leading order quantum corrections lead to the Ehrenfest 
oscillations with different power in A/a;. 

It is worth pointing out that the methods of calculating 
coupled multi-diffusons (see Sec. IV) involved in the three 
loop approximation are general. The idea may be applied 
to find the weak localization correction to other physical 
quantities (e.g. conductivity) in the ballistic quantum 
dots. 

Technically, it is convenient to use the ballistic nonlin- 
ear supcrmatrix a model (BNSoM) 30 " 32 . We will intro- 
duce this model in the next section. The quantum trans- 
port time, essential to the regularizer will be estimated. 
The supersymmctric action has two saddle points: A and 
—fcA 33 . The perturbative expansion around A gives the 
weak localization corrections to the leading perturbative 
term in the Wigner-Dyson statistics, while the pertur- 
bative expansion around —fcA gives the quantum correc- 
tions to the leading Wigner-Dyson oscillation. In Sec. Ill, 
we consider the perturbative R(w) for systems with time- 
reversal symmetry. The perturbation theory near A will 
be developed. In Sec. IV, we consider systems with bro- 
ken time-reversal symmetry. The perturbation theory 
near A will be employed to study the weak localization 
corrections in R p (co). The products of the perturbative 
expansion involve the product of multi-diffusons. We will 
calculate such coupled diffusons. Immediately, the lead- 
ing order quantum corrections give the Ehrenfest oscil- 
lations in Rp(u). In Sec. V, with the help of the global 
transformation, we will perform the perturbation expan- 
sion near —fcA. Consequently, we find the Ehrenfest os- 
cillations in the nonperturbative part of R(uj) but with 
different dependence of the period on t^. The obtained 
results are summarized in Sec. VI. Some of the calcula- 
tions of multi-diffusons are included in Appendix A. 



II. BALLISTIC NONLINEAR SUPERMATRIX a 
MODEL (BNScrM) 

Diffusive nonlinear supermatrix a model has become 
a powerful tool in the studies of R(uj) in quantum 



disorders 5 ' 33 . A natural question is where such a tech- 
nique can be generalized to the ballistic case, especially 
individual quantum dots, whether the average over im- 
purity configurations become impossible. Towards this 
direction, the first suggestion is given in Ref. 31, where 
a ballistic action was phenomenologically proposed. The 
formal justification is obtained in Ref. 30 under some cru- 
cial assumptions 34 . The complete microscopic derivation 
is given for long-ranged disorders 32 ' 35 . In this section, 
we will introduce BNScrM. Moreover, the roles played by 
the regularizer will be discussed. 



A. Ballistic supersymmetric action 

It is a standard method to introduce a supcrmatrix 
field Q. It is defined on the superspace p®g®d, where p 
stands for the advanced/retarded block, g stands for the 
fcrmion/boson block and d stands for the time-reversal 
block. Similarly to the quantum disordered case 5 , R(uj), 
Eq. (1) can be expressed as the integral over Q: 



R{w) = RcR(u), 

R(cj) = -LJdq(^J dx||STrfcAQ(x|| 



(10) 



Here the subscript in x\\ means that the integration is 
restricted on the energy shell. A is a supermatrix defined 
as: 



A : 



1 

-1 



1, 



(11) 



The ballistic supersymmctric action S is: 



s = y I ' /iSTr 



—KQ-T- 1 KCT+ — 

2 T„ 



dQ 



(12) 



where u + = w + i0 + with + being an arbitrarily small 
positive number. 1 = (r, n),rfl = drdn/2-ir. L = 
v • V — V£/(r) • ^ is the Liouville operator with U(r) 
being the classical potential. The matrix Q is generated 
byT: 



Q = T^AT, 



(13) 



and takes the value on the corresponding symmetry 
space H = G/K, where G and K are groups. For 
CUE, H = U(l,l/2)/U(l/l) <8> U(l/1). For GOE, 
H = U(2, 2/4)/U(2/2) ® U(2/2). T(l) satisfies: 



Tt(l)LT(l) = L, 
where L is a matrix defined as 



L = 



1 
k 



1, 



k 



(14) 



(15) 
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The form of the regularizer is determined by the prop- 
erties of the Liouville operator and diffraction mecha- 
nisms for particular systems. In the presence of quan- 
tum impurities with small-angle scattering, the regular- 
izer has the form presented in Eq. (12) with l/r q propor- 
tional to the impurity density 22 ' 28 . In this paper, we will 
use it to regularize the Liouville operator. According to 
the geometric theory of diffractions 36 , when an electron 
glances off a hard disk with the radius p, it emits diffrac- 
tion rays into the shadowed region deviating from the in- 

/ \ 2 / 3 

cident ray at the angle 6 <~ ( — j 37 . This leads to the 

/ \2/3 

estimation: — ~ ^f- I ) , which coincides with the 

estimation in Rcf. 38. For smooth potentials, it is esti- 
mated that Tg 1 ~ vpXp/L 2 18 - 22 , which is also obtained 
in Ref. 26 by analyzing the spreading of the wave packet. 
In any cases, r q ~ h a . It is important that the regular- 
izer is not zero but small. In Rcf. 22 (see also Sec. II D), 
it was shown that the Ehrenfest time t e — \ hi AT g . For 
\r q 3> 1, tE depends on r q logarithmically. In the semi- 
classics, r q is large but not infinite. Any change of the 
form of the regularizer accounts only for the value of the 
parameter r q under logarithm. 

B. The perturbation theory near A 

To integrate over all the modes in Eqs. (10) and (12), 
we will use the saddle point approximation. In quantum 
disorders, the Gaussian approximation leads to Altshulcr 
and Shklovskii's result 33 . For uo < t~* only the zero 



mode is important, which is coordinate independent. The 
coupling between diffusive modes accounts for the higher 
order term in the saddle point approximation. Similarly, 
in chaotic quantum dots, the the eigenmodes of the reg- 
ularized Liouville operator may contribute to R(u) in- 
dependently in the Gaussian approximation. Moreover, 
there is coupling between these eigenmodes also. 

The action, Eq. (12) has two saddle points: A and 
— fcA 33 . In Sec. Ill and IV, we will study the case with Q 
close to A. As in the disordered case, R p (uj) is dominated 
by this saddle point 30,33 . It is convenient to choose the 
parameterization below for T: 

T = 1 + iP. (16) 

The Jacobian of this transformation is unity. P anticom- 
mutes with A: 

AP + PA = 0, ^= (| o )■ (17) 

According to Eq. (14), P(l) satisfies the following condi- 
tion: 

P(l)* = -CP(1)C T , (18) 

where 1 = (r, — n) and C is the charge conjugation ma- 
trix: 

C=1^(T (19) 
The Pauli matrices: 



TO 



1 
1 



Tl 



1 

1 



Tl = 



i 



T3 = 



1 

-1 



act on the time-reversal block. Under the parameteriza- 
tion Eq. (16), we have the following expansion: 



STr [AQ] = 2^2 STr [(-l) j P 2j ] 



(20) 
(22) 



STr [JfeAQ] = 2 ^ STr [k{-l) j P 2j ] 

3=1 

s = s° + s mt - 



(21) 



and 



7TU 



J STr \p{-ioj+ + KL R )P 

/oo 
dx\\ ' 

3=1 



(-iy+ip 2 -j+i(- lL0 + + ac r )p 



(23) 
(24) 

(25) 



Here the Perron- Frobenius operator Cr is defined as 39 : 

(26) 

These expansions are essential to the perturbation calcu- 



4 = v-V-V[T(r)~_I^. 

dp T q dej) 1 



lations. 

Keeping the terms up to P 2 , we reproduce the leading 
terms of the universal Wigncr-Dyson results for R p (u>) 
[see Eqs. (27) and (28)]. In Sec. Ill and IV, we will 
calculate the weak localization correction to R p (uu). The 
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feature of finite t erg implies the existence of the gap in the 
spectrum of Perron-Frobcnius operator, which has been 



found in many chaotic systems 



,40,41 



will drop out the 



C. Quantum disorders 



Now we discuss the r q dependence of R{lo). In the 



P J+ A£rP term in the effective interaction S here presence of many quantum impurities (i.e., quantum dis 



after since it gives much smaller contribution atw<t 



erg- orders), r q becomes very small {vf/L <C 1/%). The last 
term in the action suppresses all the nonzero modes. In 
this case, one can drop out the last two terms in Eq. (12) 
and put Q to be a constant matrix 5 . In this way, the 
universal RMT results are reproduced. For GOE, 



R°(co) = 1 



A 4 



TtLU d 



— ; — - sin — 

it 4 lu 2 A duj 



~ 1 - 



1-2/, )2 



"1 


. TTLU 


,oc x 




sin — — 




UI 


A _ 


Jl t 



TTLOt 

■ sin ——at 
A 



A 4 /3 1 



2ttoA 



(27) 



For GUE, 



R u (lo) 



A 2 . 2 TTLU A 2 

» sin — = 1 - 



2tt 2 



1 — cos 



2ttuA 



(28) 



In the second line of Eq. (27), we take the limit u> S> A. Actually, for large w can be asymptotically expressed 

as 



R(lu) = R p (uj) + Rn P (uj), 



(29) 



R p {co) = 1 + 



E 



(30) 



D n +sm — ^ 



(31) 



where C n , D n and £"„ are universal numerical constants. 
R(uj) consists of two parts: One [denoted by R p (lu)] is 
perturbative in A/w and the other [denoted by R np (to)] 
is nonperturbative. It is important that expanding R p (oj) 
in terms of A/cu is a result of making the saddle point 
approximation to Eqs. (10) and (12) (see Sec. III). That 
is, we expand Q near A (i.e., T near 1) and take into 
account the fluctuations of T perturbatively. In partic- 
ular, terms like (A/u>) n ,n > 2 come from perturbative 
corrections to the Gaussian approximation in Eq. (10). 
Actually in the zero mode nonlinear a model, the action 
depends only on the parameter ui/A. As a result, the ex- 
pansion of R(u>), Eqs. (29)-(31) has no other parameter 
dependence. The reason is that in the disordered case, 
the Hikami box 23 of zero modes has dispersion only at 
lu ~ r^ 1 3> Erh, which are unimportant in the univer- 
sal region lo <^ Etk- For ui > Eth, nonzero diffusive 
modes are important and R p {lo) becomes nonuniversal 6 . 
In the disordered case, the weak localization correction 
to nonuniversal R p (uj) 6 is studied in Rcf. 11. 



D. The Ehrenfest time tE 

As in the diffusive case 5 , one may try to use the zero 
mode only and thereby drop out the free Liouville term as 
well as the regularizer in Eq. (10) in the region uj <C t~^ g 
to recover universal RMT results. However, this pro- 
cedure may miss important physics as we will present 
in this paper. To favor uniform Q, a necessary con- 
dition is that the regularizer becomes comparable to 
the Liouville term, which is of the order of ~ vp/L. 
As the initial nonuniform Q relaxes to the zero mode, 
the regularizer gets increasing ~ ^-e At due to the Lya- 

punov instability. The two terms become comparable at 
t ~ tE = \n(vFT q / L)\. Hence we conclude that only 
in the region ui <C t E , Q can be put to be a constant 
matrix in the supersymmetric action to reproduce the 
universal RMT results. 

It is important that as usual, in the region uu -C t^, , the 
diffusons become self-averaging in an individual quantum 
dot. The reason is that the time scale for the determin- 
istic chaos to become random is t erg . However, the time 
for two identical trajectories (up to the Heisenberg un- 
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certainty) to decouple (i.e., become independent of each 
other) is much longer, as we will show in Sec. V, is of 
the order of ig. Thus diffusons with the same coordi- 
nates are strongly coupled in the region i^ 1 <~ to <C t~^ g . 
The Hikami box thereby acquires an additional disper- 
sion at lu <~ t^ 1 and crossovers to the disordered limit at 
to <C tg . As we mentioned earlier, for GOE, the appear- 
ance of (A/a;) 4 term is a pure quantum effect. Indeed, 
it comes from the coupling between two zero mode diffu- 
sons. More strictly speaking, it comes from the coupling 
between a diffuson and a cooperon through Hikami box. 
In this way, it is appropriate to identify such a term as 
the weak localization correction to R p (ui). For GUE, the 
exact truncation at (A/a;) 2 term by no means implies 
the invalidity of the perturbation theory, but is due to 
the exact cancellation of the weak localization correction 
to R p (lo) arising from Hikami boxes with 4-point and 6- 
point vertex, respectively. In Sec. IV C and VC, we will 
see that the (A/a;) 4 term shows up associating with the 
Ehrenfest oscillation in chaotic quantum dots. 

In this paper, we will study R{lo) in the region A -C 
u> <C terg- For typical quantum dots, there is only one 
macroscopic scale L. The Lyapunov exponent A is of the 
same order of vf/L ~ t~^ g . Thus in the crossover from 
semiclassical to quantum picture, the Ehrenfest time ef- 
fectively opens a new region A <C t^ 1 ~ lu <C i" 1 



E. Classical limit r q — > oo 

It is important that the regularizer introduces cou- 
plings between different primitive periodic orbits. In 
the absence of diffractions or quantum impurities, r q 
goes to the infinity. In such a limit, Kogan and Efetov 
proved that Eqs. (10) and (12) are compatible with the 
Gutzwiller trace formula 32 . Furthermore one may apply 
the HOA sum rule to reproduce the leading perturba- 
tive term in R p (u>), which is of the order of 1/lo 2 [see 
Eqs. (27) and (28)]. It is important that the constant 
term 1 (so-called Weyl term) in R p (io), is pure classical 
which characterizes the phase space volume of the energy 
shell 7 . So is the ui~ 2 term, which arises from the inter- 
ference between two identical primitive periodic classical 
orbits (diagonal contribution) 15,17 . Recently, by taking 
into account the off-diagonal contributions, an oscillation 
similar to those in Eqs. (27) and (28) was found 17 . 

Thus, we conclude that in the classical limit r q — > oo, 
Eqs. (10) and (12) give the leading perturbative and non- 
perturbative terms. We obtain the same results by mak- 
ing the Gaussian approximation to Eqs. (10) and (12). 
This is not surprising. We argue that to pass from the 
Gutzwiller formula to the leading 1/lo 2 term, the only 
condition is the ergodicity, hence it is reasonable to ex- 
pect that the saddle point approximation remains good 



enough in the region A <C lo <C t~^ g . Actually, such 
a term results from the free motion of diffusons. It is 
worth noticing that so far the next to leading term ~ 
in Rp(co), as one may expect from the universal RMT 
results [see Eq. (27)], has not yet been identified in the 
framework of Gutzwiller formula. As we discussed in 
Sec. II D, it is a pure quantum effect (weak localization 
correction) . 

III. THE PERTURBATIVE R(u) IN SYSTEMS 
WITH TIME-REVERSAL SYMMETRY 

For systems with time-reversal symmetry, the weak lo- 
calization correction to R p (lu) in chaotic quantum dots 
was first studied by Aleiner and Larkin 28 . Furthermore, 
we show in this subsection that such a correction can be 
written as the second derivative of the free function at 
the two loop levels. This is important because it sug- 
gests that the result here, in principle, is possible to be 
reproduced by using the replica technique 11 ' 42 . 

A. Weak localization corrections to the perturbative 

R(u>) 

From Eqs. (10), (21) and (25), we can find that R p (lu) 
in the order of (A/a;) 3 is 




c =i(0 x = ■ = 1 
= D°(l,2) 

FIG. 1. The diagrams contribute to R p (lo) in the order of 
(A/a;) 3 in systems with time-reversal symmetry. The cross 
and the square, ~ Str [fc-P 2 ] and ~ Str \kP 4 ] , respectively, 
come from the expansion in the prefactor, see Eq. (21). The 
dot, ~ io;Str [P 4 ] , comes from the effective interaction, see 
Eq. (25). In the figure, we put A = 1. 



AR°(lj) = -^Re / DP e - s ° j dlSTr [fcP 4 (l)] f d2STr [fcP 2 (2)] - 
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-^Re J DPe- s ° J dlSTr[fcP 2 (l)] J d2STr [kP 2 {2)] • Qmti/w^ j d3STr [P 4 (3)] , (32) 

where we keep the term with j = 1 in the effective interaction S mt \ 

To calculate Eq. (32), we need to use the Wick's theorem, which implies that any even order moments of the 
Gaussian integral can be factorized as the product of the moment P 2 of the Gaussian integral. The contraction rules 
are 28 : 



2-kvP (1)M P(2) = P°(l, 2)A+STr 
and 



MA, 



+ r>°(l,2)A7STr 



MA+ 



+ D°(l, 2)A7MA+ + D°(l, 2)A+MA7 (33) 



27ri/STr[MP(l)]STr[AT'(2)] = STr (V°(l, 2)M — 2)M) Ajj~A?Ajj~ + (£>°(1, 2)M — X>°(1, 2)M) Ajj~A?Ajj~ . (34) 
Here Ajf = |(1 ± A) and M = KCM T C T K. The diffuson D°(l,2) is the solution of the regularized Liouville 



(Perron-Frobenius) equation: 



(-iw + 2) = 27r5(l - 2) 



(35) 



Z>°(1, 2; i) is the conditional probability for the particle initially at 2 to appear at 1. At the limit r q — ► oo, X>°(1, 2; t) = 
<5(Ri — R,2(i))<5(ni — 112(f)) with (R,2(£), 112(f)) the classical trajectory starting from (R,2(0), 112(0)) = (R.2,ri2). 
By using Eqs. (33), (34) and the following relation: 



we can calculate the first term in Eq. (32) as 

A x 3 

7T 

21* 



■^V°(l,3) = J d22?°(l,2)2?°(2,3), 



(36) 



dld2(P (l,2)I? (2,l)2? (l,l)) 



p u (i,i)). 



(37) 



This gives the contraction shown in Fig. 1 (1). Then using the contraction rules, Eqs. (33) and (34) we can find the 
second term consists of two parts: 



2 (^j iuj J rfld2d3(P°(l,2)2? (2,3)P (3,l)2?°(3,3)} 



(-iuj) / d3( 



d{—iuj) 



r^°(3,3) 



2?°(3,3)>, 



and 



(38) 



(-tw) / dlrf2d3(P°(3,l)2?°(l,3)I? (3,2)2? (2,3)) 



— ) (-iw) j d3( 



d 



d(—iuj) 



©°(3,3) 



d 



d{—iui) 
AR° p (u) 



D°(3,3) 



Re 



These correspond to the contractions, (2) and (3) in 
Fig. 1, respectively. Here (• • •) denotes the self-averaging 
over the energy spectrum. Adding them together, we 
find that Eq. (32) gives the weak localization correction where the free energy F^w) is 
to Rp(u)) as 



d 2 



d{—iuj) 



(39) 
(40) 



F 2 {oj) = iu J dl(V°(l, 1)D°(1, 1)). (41) 
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Actually in the r.h.s. of Eq. (41), one T>° comes from the 
diffuson and the other T>° comes from the cooperon 43 . 

P 2 (w) is the ballistic counterpart of the free energy 
in disordered systems, which is obtained by using the 
replica technique 11 . In the disordered case or in the limit 
A«w« fr^ 1 , one can justify that P 2 (w) is imaginary by 
substituting zero mode diffuson, i.e. D°(l,2) = 1/ — ito 
into Eq. (41). As a result, according to Eq. (40), AR°(u) 
vanishes. 



B. The Ehrenfest oscillations 

Now we show that for chaotic quantum dots, in the 
crossover region t E x ~ to <C t~^ g , Eqs. (40) and (41) lead 



w < tj g , 

to the Ehrenfest oscillations with the period ~ t 



E ■ 



We note that the diffusons D°(l, 1) and V°(l, 1) are 
decoupled, which means 

(V°(l, 1)V°(1, 1)) = (V°(l, 1)) (V°(l, 1)) . (42) 

Actually T>°(1,1) stands for the probability of a trajec- 
tory initiating from 1 and returning to I in the phase 
space, while T>°(1, 1) stands for the probability of a tra- 
jectory initiating from 1 and returning to 1. Such two 
trajectories start from the same place in the real space 
with opposite directions. In this way, they are distant 
parts in the phase space. Therefore, we have Eq. (42). 
In Rcf. 22, it was found that 



2?°(1,1) 



r 2 H 



(43) 



the quantum disorders. In this way, essentially it is dif- 
ferent from the Wigncr-Dyson oscillation In the low to 
limit A ~ lu <C t E , such an oscillation gradually dimin- 
ishes. Actually similar corrections to R p (lu) persist in 
higher order terms. 

Alternatively, one may Fourier transform R{tu) as the 
form factor K(t) = -J J R°(uj)e~ 1 '^ t duj before taking the 
real part. R(u) = RcR°(cu) [cf. Eq. (3)]. The RMT pre- 
dicts K{t) ~ 2t - t 2 at t <C l 5 ' 8 . We know that the 
linear term is pure classical (see Sec. HE). In fact, from 
Eq. (3) we see that K(t) = 2t- t 2 6(t - t E ) at A 2 = 0. At 
t <ts the t 2 term disappears, while at t > tE it coincides 
with that predicted by the zero dimensional nonlinear a 
model 5 . In Rcf. 44, it is proved that in the case of A 2 = 0, 
classical trajectories do not contribute to the — t 2 term 
even if the off-diagonal contributions are considered. To 
reproduce this term, the authors of Ref. 44 introduce 
some correction to the probability of self-crossing classi- 
cal trajectories. 

In this section, we prove the existence of the free en- 
ergy at the two loop level by using BNSctM [see Eqs. (40), 
(41) and (3)]. In fact, as we will see in the next section, 
at the three-loop level, the free energy does exist also. In 
this way, it is natural to expect that in general, we can 
write the perturbative part of R(lu) as the second order 
derivative of the free energy F{uu) i.e., 



R p (u>) = 



A 



Re 



d 2 



d(—iu>) : 



(45) 



As a result, 



Eq. (45) was first established for quantum disorders by 
using the replica technique 11 . We are not aware of any 
proof about it in quantum chaos. 



<P°(1,I)P°(I,1)) 



r 2 

1 2 



(44) 



From this we see that T\{lo) is the effective Hikami box. 
In the limit lu <C t E x ,r 2 (o>) — > 1. The Hikami box be- 
comes a constant, as in the diffusive disordered case 23 . In 
the crossover region, it acquires an additional dispersion 
at lo ~ t E . 

Substituting Eqs. (42) and (44) into Eq. (40), we ob- 
tain Eq. (3) 28 ' 29 . Thus R° acquires an Ehrenfest oscilla- 
tion correction with the period ~ t E x . Such an oscilla- 
tion originates from the coupling between a diffuson and 
a cooperon, which is the reminiscence of Hikami box in 



IV. THE PERTURBATIVE R(u) IN SYSTEMS 
WITH BROKEN TIME-REVERSAL SYMMETRY 

In Sec. Ill we see that the weak localization correc- 
tion to R p {io) at the two loop level originates from the 
coupling between a diffuson and a cooperon. In the bro- 
ken time-reversal systems, the cooperon modes are sup- 
pressed by the applied magnetic field. As a result, Q 
takes the components k — 0, 3 only in the d space [see 
Eq. (20)]. No longer exist the contraction rules Eqs. (33) 
and (34). Instead, we can find the contraction rules be- 
low: 



AttvP (1)M P(2) =D°(l,2)A}STr 



MA7, 



P (l,2)A7STr 



MAjj" 



+ 



2)A+r 3 STr (mA,t 3 ) + 2)A^r 3 STr (mA+t 3 ) 



(46) 



and 



27U/STr[MP(l)]STr[iVP(2)] = STr (V°(l, 2)M - X>°(1, 2)M) A 7 AA+ + (D°(l, 2)M - D°(l, 2)M) A+AA" (47) 
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by direct calculations. Then applying the Wick's theo- 
rem and the contraction rules to Eq. (32), we find that universal limit at w ^ 
it vanishes as expected. 

In the disordered case, it was first found that the lead- 
ing weak localization corrections to R p (w) in the broken 
time-reversal system are at the three loop level by us- 
ing the replica technique 11 . To our best knowledge, the 
weak localization of R p (u>) in the quantum chaotic sys- 
tems with the broken time-reversal symmetry has not yet 
been studied. In this section, we prove that for the lat- 
ter case, R p (uj) can be expressed as the second derivative 
of the free energy at the three loop level. Moreover, we 
find that the (A/w) 4 term is not exactly zero as expected 
by the universal Wigner-Dyson statistics. Instead, it is Here 



an Ehrenfest oscillation, which shows a crossover to its 



A. Weak localization corrections to the perturbative 

From Eqs. (10), (21) and (25), we can find that R p (lo) 
in the order of (A/w) 4 is 



ar$(uj) = AR; i3b (u>) + ar; m (u). 



(48) 



AR» 3b (u) = ^-Re J DPe- s ° j rflSTr [fcP 2 (l)] J d2STr [fcP 6 (2)] + 

^Rc J DPe~ s ° J dlSTr[fcP 2 (l)] j d2STr [fcP 2 (2)] • Q 
where we keep the term with j = 2 in the effective interaction S mt - . 
Ai^ >M (w) = ^Re f DPe~ s ° j dlSTr [fcP 4 (l)] j d2STr [fcP 4 (2)] + 



ITTVU) 



d3STr [P 6 (3)] 



(49) 



^Re J DPe~ s ° j dlSTr [fcP 2 (l)] J d2STr [fcP 4 (2)] 
^Re J DPe- s " J dlSTr [fcP 2 (l)] J d2STr [fcP 2 (2)] 



2 

— ilXVLO 

2 



d3STr [P 4 (3)] + 

J d3STr [P 4 (3)] J d4STr [P 4 (4)] . 

(50) 



In Eq. (50), we consider up to the second order effective interaction kept with the term j = 1. 
First, we consider Eq. (49). The first term in Eq. (49) is found to be 



AR"M 



_3 /A 
4 VtT 

3 /A 

4 V 7T 



dl( 



dld2{V°(l, 2)D°(2, 1)D°(2, 2)D°(2, 2)) 
d 



-X>°(2,2) 



P°(2,2)D°(2,2)) 



d{—iuj) 

by using Eqs. (46) and (47). This gives the contraction shown in (bl), Fig. 2. The other term in Eq. (49) gives 
Ai?;, b2 M - \ (-iw) j dld2d3(P (l,3)P (3,l)2? (2,3)P (3,2)P (3,3)) 



(51) 



3 /A 
4 



^ (-tw) J d3{ 



d 



d(—ioj) 



— £>°(3,3) 



d(—iuj) 



£>°(3,3) 



D°(3,3)), 



(52) 



and 



Afl;, M H = j ) J (lU2<i3(P"(3,l)P"(l,2)p' , (2,3)p"(3,3)P' , (3,3)) 



fi 2 



G>H^) 2 



X>°(3,3) 



r>°(3,3)£>°(3,3)}, 



(53) 



which corresponds to two ways of contraction shown in (b2) and (b3) Fig. 2. Adding them together, we find the free 
energy Fsb(^) to be [cf. Eq. (45)] 
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Fs b (u>) = ^(-M / dl{V°(l, 1)2>°(1, 1)). (54) 

Now we turn to consider Eq. (50). The first term in Eq. (50) is reduced into 

Ai£ dl (a;) = -±(*) 4 y" dld2<P (l,2)P°(l,2)P (2,l)2?°(2,l)) (55) 
by using Eqs. (46) and (47), which gives the contraction shown in (dl), Fig. 3. The second term in Eq. (50) is 

Ai?;, d2 M = 2 (J^J (-iw) J dld2d3(V°(2, 3)X>°(2, 3)D°(3, 2)£>°(3, 1)D°(1, 2)) 



- 2 (^) (-««>) y d2d3(P (2,3)P°(2,3)2?°(3,2) 



d{—iuj) 



£>°(3,2) 



(56) 



which gives the contraction shown in (d2), Fig. 3. The diagram (d3) and (d4) in Fig. 3 come from the third term in 
Eq. (50): 

AR; 43 (uj) = - (J^J (-iuj) 2 J dld2dM4(V°(3, 4)2?° (3, 4)X>°(4, 3)X>°(4, 2)X>°(2, 1)X>°(1> 3)) 



1 / A 



2 V 7T 



y d3d4(X> (3,4)£>°(3,4)£>°(4,3) 



d(—i(j) 



(57) 



2 



^ (-iuj) 2 j dld2d3d4(P (3,2)2?°(2,4)P (3,4)P (4,3)2? (4,l)P (l,3)} 
^ (-ilu) 2 J dld2d3d4(P (3,l)P (l,4)2? (4,3)I? (3,2)I? (2,4)2? (4,3)) 



1 / A 

2 



(-iuj) 2 J d3d4( 
^ (-no) 2 J d3d4( 



d 



d(—iu) 



— £>°(3,4) 



d{—iuj) 



— 2?°(3,4) 



P (3,4)X>°(4,3) 



— X>°(4,3) 



d{—iun) 



D°(3,4) 



d(—iu) 
X>°(4,3)£>°(4,3)}. 



Collecting Eqs. (55)-(58) together, we find the free energy F 3 d(w) to be [cf. Eq. (45)] 

A 2 



8^ 2 



H^) 2 / dld2(P°(l,2)P°(l,2)P°(2,l):D (2,l)). 



Hence we prove Eq. (45) at the three loop level. We 
see that Eqs. (49) and (50) give the weak localization 
correction to R p (tu), up to the order of (A/w) 4 as 



(58) 



(59) 



Re 



d 2 



d{—iuj)'- 



{F 3b {uj) + F 3d (uj)) . (60) 



B. The Couplings of Diffusons 

In this subsection we will calculate the product of dif- 
fusons appearing in Eqs. (54) and (59). 

We emphasize that it is difficult to use the standard di- 
agram technique 45 to calculate these quantities. Instead, 
it is more convenient to use the semiclassical method. 



In this section, we will follow the techniques of Rcf. 22. 
This immediately leads to the Ehrenfest oscillations in 
the perturbative R(u>). It is important that the average 
(• • •) is performed over the energy for fixed potentials, 
not impurity configurations. 



1. The couplings of two diffusons 

To calculate Eq. (59), we notice that there is the fol- 
lowing relation: 



2)2>°(1, 2)2>°(2, 1)2>°(2, 1)) 
(Z>°(1, 2)X>°(1, 2)) (X>°(2, 1)X>°(2, 1)) 



(61) 



for similar reasons as discussions for Eq. (42). Actually 
T>°(1,2) stands for the probability of a trajectory initi- 
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ating from 1 and ending at 2 in the phase space, while 
D°(2, 1) stands for the probability of a trajectory initiat- 
ing from 2 and ending at 1. Thus such two trajectories 
are distant parts and decoupled. But (D°(l, 2)X>°(1, 2)) 




FIG. 2. The diagrams with 6-point vertex contribute 
to R p (lj) in the order of (A/lu) 4 in systems with bro- 
ken time-reversal symmetry. The cross and the square, 
~ Str [fcP 2 ] and ~ Str [&P 6 ] , respectively, come from the ex- 
pansion in the prefactor, sec Eq. (21). The dot, ~ iwStr [P 6 ] , 
comes from the effective interaction, see Eq. (25). In the fig- 
ure, we put A = 1. 



can not be factorized at the region near 1 or 2, where two 
trajectories become very close to each other. 

We will study the quantity (X>°(li, 2)£>°(1 2 , 2)} instead 
of (D°(l, 2)D°(1, 2)) and put li = 1 2 in the hnal answer. 
Here we use the subscripts 1 and 2 to denote the small de- 
viations in the momentum direction and position. From 
Eq. (35), we find the motion equation to be 




• =ira X = ■ = l 

=D°a.2) 

FIG. 3. The diagrams with 4-point vertex contribute 
to R p (u>) in the order of (A/oj) 4 in systems with bro- 
ken time-reversal symmetry. The cross and the square, 
~ Str [fe-P 2 ] and ~ Str \kP 4 ] , respectively, come from the ex- 
pansion in the prefactor, see Eq. (21). The dot, ~ iwStr [P 4 ] , 
comes from the effective interaction, see Eq. (25). In the fig- 
ure, we put A = 1. 



-2iw + Ci t + Ci 2 — —JTH -4tr) 2)X>°(1 2) 2) = 2tt [J lli2 2?°(l 2 , 2) + 5i a , 2 X>°(li, 2)] . (62) 

T q 0(f> u T q OCp l2 / 



We will use (ni , Ri ) to denote the motion of the center 
of mass while ((pi,pi) to denote the relative motion 22 . 

The operator (^C± 1 + Ci^j is written as 



C c .i = v F ni 



d 



C 2 = C c 



aR 1 - w < Rl >-<*T' 



d d 2 U 

VF<Pl-^ 1 ap2 

dpi p F dR z ±A 



Pi- 



(63) 

where Rj_ denote the coordinate perpendicular to the mo- 
mentum direction. For the convenience below, we further 
introduce the change of variables: 

1/2 



z\ = In 



a\ — arctan . (64) 

Pi 



At oj <C A, the deterministic motion of the center of mass 
becomes random. Thus, the quantity £>°(li, 2)D°(1 2 , 2) 
is self-averaging over Ri,ni. In this way, we can intro- 
duce the following quantity: 

M{<f>i,p\;2,2) = (P°(li, 2)23°(1 2 , 2)) 

= /^i^°(li,2)P°(l 2 ,2). (65) 



By definition we imply that \z±\ ~ |z 2 | 1- In this 
way, we can see that the particular solution due to the 
r.h.s. of Eq. (62) is of the order of 1/Xoj and depends on 
z\ — Z2 — 0. Therefore, we can ignore the r.h.s. ofEq. (62) 
since the particular solution is much smaller than 1/u 2 
in the universal region: u<A. Then averaging Eq. (63) 
over ni,Ri and a\, we obtain 22 : 



, 9 A 2 d 2 

2w + A ^ + Y^ 



-2zi 



8 



1 



d 



2t„ dz\ V 2 dz\ 



+ 7 



M(zi;2,2) = 0. 



(66) 
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Eq. (66) exists only at the Lyapunov region, i.e., p -C 
L, <j) < l(zi < -1). At z\ ~ 0, Eq. (66) is not appli- 
cable and the solution of M does not depend on z\ any 
more. Since we are interested in the region u ~ t^, the 
coordinates 1 and 2 are independent from each other. 
Hence z\ ~ z 2 , by definition, implies that the determin- 
istic relative motion (Lyapunov asymptotic instability) 
of z\ must reach z\ = in the course of time. In this 
way, Eq. (66) must be supplemented by the boundary 
condition M(uj,Zi = 0) = M' . M' will be calculated 
below. 

The solution of Eq. (66) is 

M = w 2 (u,z 1 )M', (67) 
where w 2 (w, z-\) is found to be 22 

Ld 2 \ 2 



w 2 (u), z\) = exp 



T 



A 3 



In 



Xr q e 2zi 



■7/2 



• (68) 



up to the logarithmic accuracy. 7 < 1 is a numerical 
factor. Any other form of the regularizer with the form 



of the second order elliptic operator changes the value of 
7 only, is the Ehrcnfcst time: 



t E = ^ I In Ar 9 1 . 



In the limit z\ — > -co, M becomes 

M = r 2 (w)M'. 
At the classical limit r„ — > 00, w 2 — > as z\ 



(69) 



(70) 



But 



this is not the case for finite T q . As we see from Eq. (68), 
the classical solution W2{u, Z\) (r q = 00) thereby acquires 
a lower cutoff at z\ = —5 In Xr q . That is, we may regard 

W2(i0,— Oo) aS U> 2 (w, — 5 hi ATq). 

To find M', we repeat the procedure of deriv- 
ing Eq. (62) and study a more general quantity: 
U(w;2i,2 2 ),M' = U(w;2,2). Similarly to Eq. (66), 
?7 (05; w) is the solution of the equation: 



x a a 2 a 2 e - 2z 
2!W + A 9i + YM + ^ 



a /1 — 7 o 



dzo 



2 



+ 7 



U(Z2,L 



0, 



(71) 



where z 2 = z 2 characterizes the separation of two nearby 
trajectories centered at 2. We drop out the r.h.s. for the 
same reasons as deriving Eq. (66). To find the solution of 
Eq. (71). We repeat the procedure of Eqs. (67) to (70), 
except that the boundary condition is replaced by 



U(z2 = 0,w) 



(72) 



(V°(1,1)V°(1,1)) = (V°(1,2)V°(1,2)) = 



(75) 



We would like to mention that in the case of long- 
ranged disorders with r q = 00), an equation simi- 
lar to Eq. (66) was found by using the path integral 
technique 46 . In fact, in this case, A and A 2 were calcu- 
lated [see Appendix of Ref. 22]. It is found that A 2 ~ 4A. 



Eq. (72) means that if the deviations p\ ~ p 2 ~ L[z\ — 
z 2 = 0), then the two trajectories, connecting li and 2i, 
1 2 and 2 2 , respectively, are independent from each other. 
Thus we have (V°V '} = (V ) 2 . Similarly, we have 



M' = U(z 2 = -00)= r 2 (w)- 



1 



'(-ILU) 2 ' 

Substituting Eq. (73) into Eq. (70), we find that 

ri(w) 



(D°(1,2)D°(1,2)) = M 



(-iw) 



2 • 



At co < A, 



(73) 



(74) 



2. The couplings of three diffusons 

Now we turn to calculate Eq. (54). First, wc note 
that although the coordinates in V°(l, 1) coincide with 
each other, they are not related to each other at the re- 
gion <jj <C A, because the trajectory travels for a long 
time before returning to the same position in the phase 
space. Thus, to calculate (D°(l, 1)X>°(1, 1)P°(1, 1)) in 
Eq. (54), we will calculate a more general quantity, say 
(P°(1,2)P°(1,2)P°(1,2)). 

We proceed along the line of calculating two coupled 
diffusons. Similarly to Eq. (62), we have the following 
equation: 



-3iu> + £ij + C h2 + Ci 3 



1 d 2 



1 d 2 



1 d 2 



Tq d<j) 2 l2 T q 84>\ 



P (l!,2)I? (l 2 ,2)I? (l3,2) 



= 2tt [<5 lli2 C°(l 2 , 2)P°(1 3 , 2) + 5i 2 , 2 2?°(li, 2)P°(1 3 , 2) + S^V ^, 2)V°(1 2 , 2)] 



(76) 



Then we can separate the motion into the part of the center of mass and the relative part by introducing the change 
of variables: 



12 



r ii 
r i2 



Ri 2 — Ri, 



n, 



(77) 



Here (</>i, Ri) denotes the coordinates of the center of mass. The relative motion is described by the set of variables: 
(r\ 1 , n 2 ; 61, 62)- The transverse components of rij and ri 2 are denoted as pi and/92- Under such a change of variables, 
we rewrite Eq. (76) as 



+ c Rl + c rii + c ri2 - 



2 1 

3 To 



d 2 d 2 
de 2 + ae 2 



o 2 



D°(1 1 ,2)P (1 2 ,2)P (1 3 ,2)-0, 



t 



-vp6\ — + 



d 2 U 



dp x ' Z 1 80! ' 



3 



9p 2 ' VF dR\ x P2 'dd2 



d 



(78) 



Here we expand 6, p up to the first order. We again drop 
out the r.h.s. of Eq. (76). In Eq. (78), the form of the 
regularizer has changed compared to Eq. (62). However, 
as we discussed in Sec. II A, such a change accounts only 
for the value of r q under logarithm. Consequently, up to 
the logarithmic accuracy, t F remains unchanged. In view 
of this, we ignore the regularizer in the discussions be- 
low. To solve Eq. (78), we will proceed along the line of 
calculating two coupled diffusons. Here we will show the 
main result and give the details in Appendix A. Eq. (78) 
can be reduced as 



iiu + A 



d_ \ 2 d 2 

dx 2 dx 2 



N{x;2,2) = Q, 



H('"^+(f) 2 W s M?) 2 ). <™> 

where N(x; 2, 2) is the self-averaging of 
2?°(li,2)2? (l 2 ,2)2? (l 3 ,2) over the center of mass 
Ri,ni: 



N(x; 2, 2) = N(0, p; 2, 2) 
ee CD (li,2)2Al2,2p (l3,2)) 



cfRidni 
1 1 -no 

2irA 



(80) 



The similar equation was proposed in Ref. 26. It is 
important that here A and A2 stay exactly the same as 
those appearing in T 2 (lo) (see Appendix A). Compared 
to Eq. (66), Eq. (79) differs from it in the frequency 
(w — > |w) and the variable (z — > x). Therefore, in the 
limit T q — > 00, the homogeneous solution of Eq. (79), say 
w 3 (x;cj) is 



w 3 (x;uj) = w 2 (x; -uu) 



(81) 



supplemented by the boundary condition, which is 
w^(x — 0;u>) = 1. Up to the logarithmic accuracy, we 
obtain: 



T 3 (oj) ee lim Ws(x; v) — w 3 (— A 1 lnAr ;u') 



.'3. 9Lo 2 X 2 t E 
= exp ( -iujt E - - A2 



(82) 



Then we follow the procedure of Eqs. (67)- (74) to get: 
(P (1,2)X> (1,2)P°(1,2)> = N(x^-oo,u)- ' ^ 



(83) 



In fact, we have a general expression for n— coupled 
diffusons, say T n (ui). One is able to prove that T n (u) is 
also oscillating: ~ exp(inojtE/2) in the same spirit, which 
can be obtained in an alternative method, i.e., averaging 
the classical solution (r q — ► 00) over a minimal Gaussian 
wave packet 26 . It is in order to make the following obser- 
vations: For T n (u)), the oscillating part does not depend 
on the particular choice of the regularizer except that t e 
depends on r q parametrically. 



C. The Ehrenfest oscillations in Rp(co) 



F 3b (oj) and F 3d (u) are the ballistic counterparts of the 
zero mode approximation of the free energy in disordered 
systems. In the disordered limit or A <C to <C t E , one 
can verify that Fsb(u>) and Fsd{w) cancel each other and 
Eq. (60) gives vanishing result for ARp(uj) by putting 
T>°(1, 2) to be zero mode diffuson, i.e. T>°(1, 2) = 1/ — iuj. 
However, as we will see below, it is not so for chaotic 
quantum dots in the crossover region t^} ~ u -C t~^ g . 
Instead, Eqs. (54), (59) and (60) lead to the Ehrenfest 
oscillation with the period ~ t E x . 
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Substituting Eq. (83) into Eq. (54), we find that 

A 2 rlM 



8tt 2 (-iw) 2 ' 



(84) 



Taking Eq. (61) into account, we find 



F 3d {co) 



a 2 r|H 

8tt 2 (-zw) 2 ' 



(85) 



Therefore, we obtain Eq. (7). We see that -Rp(w) ac- 
quires an oscillation correction at Lots ~ 1 with the pe- 
riod ~ i^ 1 . 

We should emphasize that the appearance of the 
Ehrenfest oscillations, ~ e mujtE ,n = 1,2,3,- •• do not 
depend on a particular form of the rcgularizcr, although 
tE does. This can be understood in the following way. 
In the Lyapunov region, n coupled diffusons (Hikami box 
with 2n-point vertex) contribute to R(iv) (in the disor- 

dered limit) an additional factor ~ ( Gf +i!L G^_^_ J . In 

the semiclassics, Gf ~ e ul ' VF and Gf ~ e-" J /" p with I 
being the length of the trajectory piece. Up to the log- 
arithmic accuracy, the Lyapunov regions have the same 

length v F t E . In this way, (g*_„ G?_„ )" - e in " tE . In 

general, R p (ui) can be expressed as 

i?» = l + W-) C n {ut E ) 



= 1 + 



n 

a 2 d 2 

it 2 doj 2 



F(w) 



A 



L-l 



E^+y-s,iCsMB), (86) 



where L, V and 5 characterize the topology of the di- 
agram of the free energy, from which R p (u>) is ob- 
tained. L, S and V are the number of loops, sides (dif- 
fusons) and vertices (Hikami boxes) respectively. The 
constraint comes from the well known Euler theorem in 
the topology 47 . For a topological diagram of the free en- 
ergy with S sides (diffusons), we have 
e iSLot E j n j- ne systems with time- reversal symmetry, the 
leading corrections come from the two loop diagram, 
L = 2,S = 2 and V = 1, hence the leading Ehrenfest 
oscillation in R p (u>) is ~ e l2uJtE . In systems with bro- 
ken time-reversal symmetry, the leading corrections come 



from three loop diagrams. The diagram with the Hikami 
box of 6-point vertex (L = 3, S = 3 and V = 1) leads 
to the Ehrenfest oscillation ~ e i3a)tE . The diagram with 



two Hikami boxes (L 
Ehrenfest oscillation ' 



3,5 = 



4 and V = 2) leads to the 



V. THE NONPERTURBATIVE PART OF R(u) IN 
SYSTEMS WITH BROKEN TIME-REVERSAL 
SYMMETRY 

In Sec. Ill and IV, we study the perturbative expan- 
sion around the saddle point A. As a result, it gives the 
weak localization corrections to the leading perturbative 
term of the universal Wigncr-Dyson statistics. However, 
it is known that such an expansion cannot reproduce the 
Wigner-Dyson oscillation 33 . Instead, Andreev and Alt- 
shuler found that such a nonperturbative term is con- 
trolled by other nontrivial saddle points 33 . In this sec- 
tion, we apply their method to study the leading order 
corrections to the universal Wigner-Dyson oscillation in 
the region A«u~ t E x <C t~} g . 



A. The global transformation 

For GOE and GUE, the other saddle point is -fcA 33 . 
To establish a perturbation theory, a trick, so-called 
global transformation on H is introduced, which maps 
the saddle point — kA to A. That is, we perform a global 
coordinate transformation: 



r = U^QUo, T^T' = TC/ , 



(87) 



where Uq € H satisfies — Uq kAUo = A. It is important 
that the Jabobian is unity and Uq does not depend on 
the coordinates. Under the transformation Eq. (87), 



Q' = T'^AT'. 
Thus, we can rewrite R(oj), Eq. (10) as 



(88) 



R(u) = — Rc I DQ'[ I dx,|STr [AQ'(x,|)] ) e 



where S e ff is the action: 



(89) 



iris f 



(foiiSTr 



0(p 



(90) 



Although the findings in Rcf. 33 are asymptotic, it 
is important that the perturbation near the nontrivial 
saddle points take into account nonzero mode approx- 
imations, which, in principle, cannot be done by zero- 



dimensional supermatrix a model. Indeed, for BNSctM, 
the Gaussian approximation is not applicable 32 in the 
nonuniversal region: ui 3> t~^ g . However, we believe 
that the perturbation near saddle points is applicable 
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in the universal region u> <C t e ^ g , because the repetition 
problem 17 is not essential in this case. In the ballistic 
case, the new region: t^ 1 ~ uj <C t~} g , appears, within 
which the relaxation of the momentum is not complete. 
As a result, one can not put Q to be a constant matrix. 
The parameterization, Eqs. (16)-(18) still holds for T". 
Furthermore, we can expand B in terms of the Pauli's 
matrices: 

B = ^-B fc ®r dife , B = ^B k ®T d . k . (91) 

k k 

The modes k = 0, 3 stand for the diffuson and the modes 
k = 1,2 stand for the cooperon. They satisfy the same 
equation Eq. (35) in the absence of magnetic fields. In 
this way, we call both of them diffuson modes in this pa- 
per unless special explanation. The modes k = 1,2 are 
suppressed in systems with broken time-reversal symme- 
try due to the destruction of the interference between a 



trajectory and its time-reversed partner. In this paper, 
we will study only this case. Due to Eq. (14), the matri- 
ces Bk and Bk satisfy: 

B k = kBl fc = 0,l,2,3. (92) 

Bk may be parameterized as 

*-(<*)• <93 » 

where afc(l) and fefc(l) are ordinary variables, while Ofc(l) 
and rjk(l) are Grassmann variables. Then wc substitute 
Eqs. (91)-(93) into Eq. (89) and keep the quadratic terms 
in cifc, bk, (7fc and i)* k to rewrite Eq. (90) as 

S° eff = -2i7TLu + S°. (94) 



S° = 2ir dx\\ [ a U x \\) [ iuj+ +£r) ak( x \\) + K( X \\) (-«^ + + £r) h(x\\) + aU x \\)£R°k(x\\) + m(x\\)£ R r)U x \\) 

J fe=0,3 

(95) 



It is important to see that the new Gaussian action, S° 
does not include zero mode contribution of Grassmann 
fields. Thus, for a generic Gaussian integral with the ac- 
tion S° not to vanish, the prefactor must consist of all 
the independent zero mode Grassmann variables. In this 
way, by expanding Q' in the prefactor up to the second 
order, and taking into account the zero mode bosonic 
and Grassmann fields only, the universal Wigner-Dyson 
oscillation is reproduced in Refs. 30, 33. 



X i X I 




(2) 




(3) 



X i X 



: ICO 



= ■ = 1 



B. Quantum corrections 



D°(l,2;o>) 



= D°(l,2;0) 



The leading order corrections to the universal Wigner- 
Dyson oscillation has the amplitude ~ (A/ui) 3 at uj 3> A. 
To see this, one need to expand T" -1 and keep Q up to P 4 
and take into account the first order effective interaction 
similarly to Eq. (25). As a result, we obtain: 



FIG. 4. Nonvanishing diagrams contribute to R np (ui) in the 
order of (A/a;) 3 for systems with broken time-reversal sym- 
metry. The cross and the square, ~ Str [P 2 ] and ~ Str [P 4 ] , 
respectively, come from the expansion in the prefactor, see 
Eq. (22). The dot, ~ iojStr [fcP 4 ] , comes from the effective 
interaction. The solid line stands for the bosonic propagator, 
the dashed line stands for the zero mode Grassmann propa- 
gator, and the dash-doted line stands for the nonzero mode 
Grassmann propagator. In the figure, we put A = 1. 



8 / rflSTr 



$ifj = -^-iuJ I dx||STr[fcP 4 (:r|| 



[P 2 (l)] J rf2STr [P 4 (2)] - 4 J dlSTr [P 2 (l)] J d2STr [P 2 (2)] S l e f f 

(96) 
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To calculate Eq. (96), we need to separate the zero mode Grassmann fields from all the other fields and integrate 
them out first. Consequently we are left with a Gaussian integral involving the bosonic and nonzero mode Grassmann 
fields. Then we can factorize the integral into the product of the following pairs: 

1-kvZ- 1 j D[a*ab*bo*ariri*]a* k {l)a k {2)e- §0 = 2; -u), (97) 

AnvZ- 1 J D[a*ab*b<j*<j m *]b* k {l)b k {2)e-~ sa = V°(l, 2; w), (98) 

A-kvZ' 1 j D[a*ab*ba*ar]r]*]a* k (l)u k {2)e-' s0 = X>°(1, 2; 0), (99) 

47TI/Z- 1 J L>[a*a6*6crV W *]r/ fc (l)77^(2)e- s0 = P°(l, 2; 0), (100) 



where k = 0,3, and 2?°(1,2;0) satisfies the following obtain Eq. (102). The measure D [a* ab*ba* or\rf\ is on all 

equation: the independent bosonic and nonzero mode Grassmann 

fields. Besides, one needs to take into account Eq. (14) 

£rV (1, 2; 0) = 2n [5(1 - 2) - 1] . (101) to reduce the number of integration variables by half. 

In the square integrable space, the overall factor, Z can Eq. (96) involves the Gaussian integral of bosonic and 

be approximated as nonzero mode Grassmann fields, which can reduced into 

j the sum with each term the products of pairs [Eq. (97) to 

D[a*ab*ba*a7]7]*]e~ s ° ~ (102) (100)]. Eq. (96) corresponds to 4 nonvanishing diagrams 

- (4tt)2 ^ 

shown in Fig. 4. One should not confuse these diagrams 

„ T ,.1.^/11 an ■ with Fig. 1. Here the Hikami box in the 8-shaped dia- 

Wc note that Z =t= 1 because b is not supersymmtric , , . ,, .. . ,.rr 

, , . ' , „ „ , , gram stands tor the coupling between two diliusons, not 

and we do not integrate out zero mode Grassmann fields. , . , TJ 

a ditiuson and a cooperon, as m the orthogonal case. It 



Actually with nonzero modes integrated out, we come 
up with a factor (1 + a), where a = o-i 2 /^2 k A| and the 



is important that such a diagram is exactly zero in the 
unitary case without the global transformation. In this 



sum is over all the nonzero eigenvalues X k of the Perron- ,. . , ,. ., „ ,, 

^ . lrnnls °„ . . . way, the existence ol these diagrams results irom the su 

Irobcmus operator ' ' In the universal region uo <C , , . 

A, this factor crossovers to the universal value 1. Further- 

more, after integrating out zero mode bosonic fields, we The first term in the r.h.s. of Eq. (96) gives 



Ai£ pil (w) = (e 2 ™/ A J d3 [{D°(3, 3; -w)) - (X>°(3,3;w)>]^ ~ -4-^Ime 2 " w / A (103) 

which is the contraction shown in Fig. 4(1). Above we use the fact that at w < A, we need to take into account the 
zero mode only in the diffusons. 

The second term in Eq. (96) has two nonvanishing contraction. One is shown in Fig. 4(2), while the other is shown 
in Fig. 4(3) and (4). In Fig. 4(2), only the bosonic fields are involved. By using Eqs. (97) and (98), we find it to be: 

Ai?^ 2 M = -^Rc (V^) e 2 ™" / A | d2 j d3[(P (3,2;-c)2? (2,3;-c)) + (P (3,2; W )P (2,3; W ))]) . (104) 

It should be emphasized that the order of the arguments appearing in the product of two diffusons above is essentially 
different from what we calculated in Sec. IV B 1 [see Eq. (74)]. In fact, we can use Eq. (36) to reduce the two diffusons 
in Eq. (104) into one diffuson. Thus, we can approximate the product of two diffusons as l/(±iu)) 2 in the region 

Fig. 4(3) differs from Fig. 4(4) in that there the coupling occurs between bosonic fields. While Fig. 4(4), the coupling 
occurs between the bosonic and nonzero mode Grassmann fields. By using Eqs. (97) and (98), we find that Fig. 4(3) 
corresponds to 

Ai£ pi3 (a;) = -^Rc (Vzc)e 2 — / A J d3 [(P°(3, 3; -u)V°(3, 3; - W )> + <2?°(3, 3; lo)V°(3, 3; W )>] ) . (105) 
By using Eqs. (97)-(100), we find that Fig. 4(4) corresponds to 

Ai?^ 4 H = ^Rc ({-iu)e 2 ^ ' A J d3 [<2?°(3,3;- W )2?°(3,3;0)) + (V°(3, 3; lj)V°(3, 3; 0))]) . (106) 



16 



In the disordered case, one is easy to check 
that AR^ p (u) vanishes using the zero mode diffuson: 
(X»°(l,2;w)) = 1/ -iuj+ and (X»°(l,2;0)) = 0. However, 
it is not so in the crossover region A«w~ t^ 1 -C t~^ g . 
Instead, as we see below, Eqs. (103)-(106) lead to the 
Ehrenfest oscillation with the period (A -1 + atE)~ 1 with 
a a universal numerical factor. 



C. The Ehrenfest oscillation in R™ p (lo) 

Now we show that in the crossover region, Eqs. (103)- 
(106) lead to the Ehrenfest oscillation, but with different 
dependence of the period on tE compared to the pertur- 
bative part. 

First, according to Eq. (75), we can rewrite Eq. (105) 

as 



A 3 



-Re 



-iuj 



(107) 



In the case of A 2 w 2 /A 3 <c 1, Eq. (107) is simplified as 

AR" 3 (w) = -2^sin^cos2u;t E . (108) 

To calculate Eq. (106), first we under- 
stand £>°(3,3;0) in the way of £>°(3,3;0) = 
lim wi ^ [£>°(3, 3; w{) + X>°(3, 3; -wi)] /2, where 
P°(l,2;wi) should be considered as the solution of the 
equation below: 

(-iwf + C R ) 2; wi) = 2tt [S(l - 2) - 1] . (109) 



Then following the procedure of deriving Eq. (74), we 
find that (D°(3, 3; w)£>°(3, 3; wi)) = -r 2 
Thus we obtain: 

(P°(3,3;c)P°(3,3;0)) = -^r* (f ) • (HO) 

Substituting Eq. (110) into Eq. (106), we find that 
Ai£ Pi4 (w) is 

A 3 27TO; 9 n LlJ 

ARl pA {u) = -2— sin— ^Rc— r 2 (-). (Ill) 
In the case of A2^ 2 /A 3 -C 1, Eq. (Ill) is simplified as 



A 3 27TW 

2 sin — — ujtE sinwig. (H2) 

7TUr A 



Collecting Eqs. (103), (104), (107) and (111) to- 
gether, we find that the leading quantum corrections in 
the nonperturbative part of R(u>) is Eq. (9). Thereby 
Rn P (uj) acquires an oscillation correction with the period 
(A -1 + atE)^ 1 with a a universal numerical factor. Ac- 
tually we can find the similar results for systems with 
time-reversal symmetry, but the amplitudes are much 
smaller, which are proportional to (A/ui) 5 cos(27ru>/A) 
or (A/w) 5 sin(27ra>/A). 

In the low u limit: ujt E < 1, AR^ p (uj) ~ A 3 i|/w. It 
is a small quantum correction to the universal Wigner- 
Dyson statistics. However, we are not able to find this 
term by using the method developed in Ref. 10. In 
the large to limit: u)tE 3> 1, we see that the universal 
Wigner-Dyson statistics acquires a small universal cor- 
rection ~ (A/w) 3 e 27ri "/ A . But it is smaller than the 
nonuniversal correction, which is ~ (A/A) 2 e 27rlW//A . 

Taking into account Eqs. (86) and (9), we can formally 
express R{uj) as the following : 



C0S ^-Z^( — ) D n(ut E ) +sm — 2^ ( — j E n (wt E ) 



(113) 



The factor a ~ w 2 /A 2 arises from making the Gaus- 
sian approximation to the nonzero modes of the Perron- 
Frobenius operator 30 ' 48 . We do not study this term 
in this paper because it does not lead to the Ehren- 
fest oscillations. For disordered systems, a similar term, 
~ uu 2 /E^h was found in Rcf. 10. In general, C n ,D n and 
E n are oscillating functions of wis and proportional to 
T n (u>). In the disordered limit, C n ,D n and E n become 
universal constants. 



VI. CONCLUSIONS 

In this paper, we study the electron energy statistics in 
chaotic quantum dots with one macroscopic size L. For 
such systems, the inverse ergodic time and the Lyapunov 



exponent are of the same order: t~^ g ~ A ~ ^£-. Con- 
sequently, in the semiclassical limit, the Ehrenfest time 
tE opens an intermediate region A <C t E l <~ uj <C t~^ g . 
We study the behavior of the two level correlation func- 
tion, R{uj) in the universal region A«wC t~} g in sys- 
tems with broken time-reversal symmetry. Surprisingly, 
it is found that R(u>) deviates from the universal Wigner- 
Dyson statistics. Basically, we find that R(u>) acquires 
two types of oscillation corrections in the crossover re- 
gion t E x ~ ui <C t~r g - The oscillation periods have differ- 
ent dependence on tg. The Ehrenfest oscillations in the 
perturbative part have the period proportional to t^ 1 , 
while the Ehrenfest oscillations in the nonperturbative 
part have the period (A -1 + ats)^ 1 with a being some 
universal numerical factor. These additional Ehrenfest 
oscillations are small corrections to the universal RMT 
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results. In particular, for Rp(ui) described by Eq. (28), 
the exact truncation at the term ~ oj~ 2 does not imply 
the disappearance of quantum corrections. Instead, it 
is due to the cancellations of contributions arising from 
Hikami boxes associating with different kinds of vertices 
at the limit uitE — > 0. The Ehrenfest oscillations are just 
the reminiscence of such cancellations. 

In this paper, only the leading Ehrenfest oscillation 
corrections are calculated. Actually there is a general 
expression, Eq. (113) in the region A <C u> <C t~^ g . The 

functions C n (x) [hence C n {x)\, D n {x) and E n (x) are os- 
cillating at u>tE ^ 1- In the limit nojtE C 1, C n (hence 
C n ), D n and E n crossover to their universal values. In 
particular, in the GUE case, C 2 = — D 2 — —1/2, and all 
other coefficients vanish. The universal values of C n {x) 
[hence C n (x)], D n (x) and E n (x) do not depend on the 
magnitude of the regularizer [the last term in Eq. (12)]. 
But it depends on the order of how we take the limits: 
cj — > and r q — > 00. Actually we take the limit ui — > 
before r q — > 00. In contrast, the boundary of the region 
where the universality exists [u <C t^ 1 ) does logarith- 
mically depend on the magnitude of r q . The regularizer 
describe the coupling between classical trajectories. In 
this sense, it plays a role similar to that played by the in- 
teraction between particles in nonidcal Bose gases. There 
the scaling theory describes the universal behavior at the 
transition point and the critical indices do not depend on 
interactions. However, the region for the scaling to be ap- 
plicable does depend on the interaction. While in ideal 
gases, the scaling is absent. It is important that the re- 
sults presented here hold only at A -C oj <C A, where 
the saddle point approximation to Eqs. (10) and (12) is 
possible. 

We note that T n (ui) is not analytical at l/r q in the 
large uj limit tg 1 <C u) < A. Actually following from 
Eqs. (4) and (82), we find that 



r„(w) ~ — 



u 2 /(\ 3 /\ 2 )+iw/\ 



t E x < w < A. (114) 



t ej r g , the nonzero mode contribu- 



From this, we see that corrections to R(oj), predicted 
by the Gutzwiller formula, are proportional to T~ a ,0 < 
a«l. In this way, any attempts of establishing a l/r q 
expansion are prohibited. 

At higher uj, uj > A 
tions turn out to be important. In this case, the saddle 
point approximation is no longer applicable. A refined 
technique is desired to study the oscillation, claimed in 
Ref. 48 for generic chaotic quantum dots. 

In this paper, we consider chaotic quantum dots where 
A ~ tern ~ ^ • Actually, it is possible that A 3> . In 



the latter case, the dot has large enough size to contain 
a lot of classical impurities inside. In this case, the Lya- 
punov exponent A ~ navp, where n is the concentration 
of impurities and a is the size of classical impurities 49 . 
The ergodic time (Thouless time) t erg = L 2 /D 3> A -1 . 
Here D is the classical diffusion constant. In this case, 



contributions from nonzero modes are possible to be de- 
scribed by the saddle point approximation. The similar 
Ehrenfest oscillations are expected to exist at t~^ g < uj ~ 
i^ 1 < A. We leave this work in the future. 

Based on the present work, it is unclear whether the 
BGS conjecture 13 may work in the universal region uj <C 
terg f° r generic quantum chaotic systems. In other words, 
are the fluctuations of energy levels described by the uni- 
versal Wigner-Dyson statistics in such systems? Accord- 
ing to Eq. (113), we point out that in chaotic quantum 
billiards a new scale, the Ehrenfest time appears, which 
results in the new types of oscillations at to ~ t E . These 
oscillations arc not expected by either the RMT, zero 
mode nonlinear supermatrix a model or the Gutzwiller 
trace formula. Secondly, despite of the appearance of the 
Ehrenfest oscillations in the region uj ~ t E x 3> A, their 
amplitudes are small. In this way, R(uj) is still domi- 
nated by the universal Wigner-Dyson statistics. Finally, 
the behavior of R(lu) at uj ~ A remains an open prob- 
lem. Indeed, in the present work we use the saddle point 
approximation to take into account the zero mode con- 
tributions. Furthermore, the perturbation theory near 
the saddle points predicts the regions for the terms in 
the expansion Eq. (113) to approach their universal lim- 
its. However, we are not aware of the limit of the sum of 
such an asymptotic expansion. 
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APPENDIX A: THE DERIVATION OF 
EQUATION (79) 

In this Appendix we prove Eq. (79) . We will follow the 
general method developed in Ref. 22. In the discussions 
below, we ignore the small regularizer in Eq. (78) for the 
moment. Moreover, we turn to the time representation. 
Then employing the change of variables 



z = (21,22), a = (ai,a 2 ), 



Z\ = In 
Z2 = hi 



? 2 2 + 



Pi 



(El 



1/2 
1/2 



a.\ = arctan 



1 

Pi 



0L2 = arctan , (Al) 

Pi 



we rewrite Eq. (78) as (after averaging over the coordi- 
nates of the center of mass) 



d_ 

di 



+ D 1 +D 2 ) N(z,a;2,2) = 0, 
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b x = -Bi(t)sin2ai— + [Bi(t) cos2oi + B 2 (t)] — , 
az\ aa.\ 

d d 
D 2 = -B x (t) sin 2o 2 — + [Bi (t) cos 2o 2 + B 2 (t)} — , 



parametrically and has no ot\ %2 dependence. The evo- 
lution of N(t; x, y; 2, 2) is governed by the Fokker-Planck 
type equation: 



VF 



T 77 



' dz 2 
L d 2 U 



2L ' 2p F <9i? : [ |R - R(t ' Ro) ■ 
Introducing the change of variables 

zi + z 2 



x = 



y = z x - z 2 , 



(A2) 



(A3) 



dt 



d_ 

\dx 



N(t;x,y;2,2) = 0, 



where T (^) is defined as 



T 



(- 

\dx 



= lim - In 

t— »oo t 



dndR 

2nA 



exp 



B(t) 



d_ 

dx 



(All) 



(A12) 



we obtain: 

Di+D 2 = 
1 



Here we keep in mind the important fact that 

J* dtBi(t) sin [2ai y2 (af 2 , t)] does not depend on oi j2 
gdue to the self- averaging at large times, so that 



--Bi(t) (sin2oi +sin2o 2 ) — + B^t) (sin 2a 1 - sin 2a 2 ) —r* dtBl sin [2di(a?,t)] = /* dfcBi (t) sin [2a 2 (a§,t)]. 



+ [Si(t) cos 2 ai + B 2 (t)] — + [B x (t) cos 2a 2 + B 2 (t)} - 

ooii oa 2 



q In this way, T has no dependence. 



The formal solution of Eq. (A2) is 
N(t;x,y;2,2) = 



exp 
where C\ is 
Ci (t,ai, 2 ) 



Ci {t,ai j2 ) 



d_ 

dx 



N(0;x,y,a°(a,t);2,2), (A5) 



(A4) B(t)= [ dtBi (t) sin [2&(a°,t)] , (A13) 

Jo 

where a can take the subscript cither 1 or 2. 

Because N(t; x, y; 2, 2) slowly varies over x, we expand 
T up to the second order in ^ to get: 



T 



(d\ ^ x d_ + ^_&_ 

\ dx J dx 2 dx 2 ' 



2io 



dtiBi(ii) (sin2di + sin2d 2 ) . (A6) 



The functions of y(t) and 01,2(0° 2 ,t) satisfy the follow- 
ing equations: 



where 



and 



lim 

t— >oo 



7/ 



2ttA 



B(t) 



<9y 

at 



= Bi(t) (sin2oi — sin2o 2 ) , 



d 

Wt 



h,2 = Bi(t)cos2di, 2 + B 2 (t), 



"1,2(01,2,0) = a 1>2 , 
with di j2 implicitly defined by 

"1,2 [«l,2( a l,2)*)i*] = Q!l,2- 

The formal solution of y(f) is 



(A7) 

(A8) 
(A9) 



lim -{ 

t^oo t 



I 



dndU 

2ttA 



B 2 (t) 



- AT}. 



(A14) 



(A15) 



(A16) 



It is important that here A and A 2 are exactly the same 
as those appearing in T 2 (co) because B(t) is the same as 
that in Ref. 22. This is not surprising because mathemat- 
ically, A is determined by the eigenvalue of the stability 
matrix of a trajectory 50 . Returning to the frequency rep- 
resentation, we find Eq. (79). 



y{t)=y(0)+[ dtiBi(ti) (sin2oi - sin2o 2 ) (A10) 
Jo 

with 2/(0) being the initial condition. 

We are interested in the region t > A -1 . In this case, 
di, 2 (di 2 ,t) become self-averaging over the coordinates 
of the center of mass and no longer depend on the ini- 
tial conditions. Consequently, at such large times, y(t) 
reaches a constant depending on the initial value [see 
Eq. (A10)]. In this way, N(t;x,y;2,2) depends on y 
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